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EFFECTS OF STRUCTURAL DAMPING ON 
FLUTTER OF STRESSED PANELS* 

By Charles P. Shore 
Langley Research Center 

SUMMARY 

The flutter of stressed orthotropic panels with edge rotational restraint is investi- 
gated theoretically. A modal solution which utilizes two-dimensional quasi-steady aero- 
dynamic theory and includes both structural and aerodynamic damping is presented. The 
investigation was conducted in an attempt to explain the existing discrepancies between 
theoretical and experimental results for stressed panels, in particular for panels stressed 
near buckling, the most critical flutter condition. Structural damping is represented in 
the present analysis in a manner consistent with the representation for a Kelvin-Voigt 
viscoelastic body. Numerical results are presented to indicate that such a representation 
of structural damping eliminates the physically untenable results that have been obtained 
for certain values of the input parameters in previous analyses. Additionally, represen- 
tation of structural damping in this manner results in a reasonable correlation between 
theory and experiment over the entire range of stress from zero to buckling. 

INTRODUCTION 

Recent studies of panel flutter (refs. 1 and 2) have shown that inclusion of the effects 
of edge rotational restraint can improve the agreement between theoretical and experi- 
mental flutter results for stressed isotropic panels. However, many instances remain 
where large discrepancies exist between theory and experiment for such panels. Theory 
which neglects damping predicts zero dynamic pressure required for flutter whenever the 
midplane load has caused two panel vibration frequencies to be equal. This condition 
occurs for panels with a length-width ratio greater than 1 for many combinations of stress 
ratio and edge rotational restraint. On the other hand, experimental flutter boundaries 
display none of the anomalous behavior associated with coincidence of the theoretical 
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frequencies. In fact, most experimental data show that the minimum value of dynamic 
pressure for flutter occurs at or near the transition point, which is the intersection of the 
panel prebuckled boundary and the postbuckled boundary. For examples of this behavior, 
see references 2 and 3. 

Inclusion of aerodynamic and structural damping in the flutter theory can remove 
zero-dynamic-pressure flutter points. (See refs. 4 and 5.) However, the effects of aero- 
dynamic and viscous structural damping (employed in refs. 6 and 7) decrease with 
increasing compressive stress and vanish at the transition point (ref. 5). Reference 8 
points out that viscous-type damping is not a realistic approximation of the damping char- 
acteristics of structural materials and should, therefore, be applied with caution. Another 
formulation for structural damping frequently utilized in aeroelastic studies is the fre- 
quency independent linear hysteretic form. Although the mechanisms responsible for 
damping in structures are usually nonlinear, for small amounts of damping the linear 
model often provides a good approximation (ref. 9). In references 5 and 10 to 14 a linear 
hysteretic structural damping coefficient is, in effect, used to modify both the bending and 
membrane loading terms of the flutter equation. However, the effect of structural damping 
when employed in this manner also decreases with an increase in compressive stress and 
vanishes at the transition point (refs. 4 and 5). Hence, this approach predicts the same 
anomalous behavior as the theory that ignores damping. However, vibration equations 
derived on the basis of viscoelastic theory utilizing stress -strain relations for a Kelvin- 
Voigt body (ref. 15) indicate that linear hysteretic structural damping should modify only 
those terms of the equations associated with bending. Such a procedure is believed to 
lead to the correct form of the equations and is the basis for the flutter analysis developed 
in the present paper. However, in order to compare present results with results from 
previous investigations, two separate damping coefficients are introduced, one associated 
with the membrane or inplane forces and the other associated with panel bending. A solu- 
tion for the flutter of orthotropic, flat, rectangular panels with arbitrary edge rotational 
restraint is obtained by application of the Galerkin procedure. It is assumed that the 
lateral loading is given by two-dimensional quasi-steady approximate aerodynamic theory. 

Results from the analysis are presented to indicate that the use of the two- 
dimensional quasi-steady aerodynamic theory yields satisfactory results for problems in 
which the effects of aerodynamic damping are important. Additional results are presented 
to indicate that inclusion of structural damping in flutter theory is more realistically 
accomplished by use of linear hysteretic structural damping associated with panel bending 
only. Such a procedure removes all physically untenable effects that otherwise occur and 
renders the value of the flutter parameter at the transition point relatively insensitive to 
variations in length-width ratio, stress ratio, and edge rotational restraint. Finally, a 
comparison of results from the present analysis with existing experimental results is 
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presented to verify the improved correlation between theory and experiment for flutter of 
stressed panels over the entire stress range from zero to buckling. 

SYMBOLS 


A rs Fourier series coefficient 

a panel length 


Bj., B 2 , B 3 , B 4 coefficients defined by equations (17) 


b 


C 0 , ^2 

c 


n 


X 


D t 


D xy 

D 1 > D 12> D 2 

E 


F r 

Gs 

g 

g A 

§B 


panel width 

constants defined by equations (A2) 
free-stream speed of sound 
panel bending stiffness in x-direction 
panel bending stiffness in y-direction 
panel twisting stiffness 

panel stiffness coefficients defined by equations (9) 
Young’s modulus 
streamwise deflection function 
cross -stream deflection function 
structural damping coefficient 

pc 

aerodynamic damping coefficient, 

yu r 

bending structural damping coefficient 


i = fl 


membrane structural damping coefficient 
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(Il)mr> (I 2 ) r ’(*3) mr coefficients defined by equations (A3), (A4), 

( J 2 ) mm -( J 2 )mr’( J3 ) mm coefficients defined by equations 

(Kj) ,(Ko) , (Ko) , (K 3 ) coefficients defined by equations 

' 'mm ' 7 mm ' /mr ' T 


mm 


K x,cr 

k x ,T 

k y 

ki,k2 

^mr 

k'mr 

Vr 


M 

N x 

Ny 

p 

q 


N*b 2 

nondimensional stress coefficient in x-direction, — — 

7t2Di 

buckling load with no airflow 
buckling load with airflow 

nondimensional stress coefficient in y-direction, 
constants 

coefficient defined by equation (A 6 ) 
coefficient defined by equation (A20) 
coefficient defined by equation (A21) 


Nyb 2 

7t2Dj 


l (x,y,t) aerodynamic loading term 


Mach number 

inplane loading in x-direction, positive in compression 
inplane loading in y-direction, positive in compression 
modulus of a complex number defined by equation ( 22 ) 
free-stream dynamic pressure 

rotational restraint coefficient on x = ±§- boundary, — — 

2 Dj 

rotational restraint coefficient on y = ±j 2 boundary, — ^ 

2 D 2 


and (A 5) 

(A12) to (A15) 
(A16) to (A19) 
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t time 

U free- stream air velocity 

w lateral deflection of panel 

x,y,z Cartesian coordinates of panel 

a. complex frequency coefficient, cp + iu; 

/3 compressibility factor, \Jm^ - 1 

y panel mass per unit area 

6j constant dependent on degree of rotational restraint 

e strain 

€q maximum strain amplitude 

6 x ,9y spring constant of rotational spring supporting panel on boundary x = ± 

and boundary y = ±j|, respectively 

X dynamic-pressure parameter 

dynamic-pressure parameter at transition point 

M x ,^y Poisson T s ratio in x-direction and y-direction, respectively 

p free-stream air density 

a stress 

cp real part of complex frequency exponential coefficient 

i p complex eigenvalue (see eq. 18) 

\Jsjl real part of complex eigenvalue 

imaginary part of complex eigenvalue 
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to IP 



CO 


circular frequency 


u> r reference frequency 

Subscripts: 

m,r number of half-waves in streamwise direction 

s number of half-waves in cross-stream direction 

ANALYSIS 


The flat, rectangular, orthotropic panel configuration under consideration is shown 
in figure 1. The panel which is of length a and width b is exposed to supersonic flow 
over one surface and is subjected to uniform inplane force intensities N x and Ny which 
are positive in compression. The panel is supported such that the lateral deflection on 
all edges is zero. Additionally, the panel edges are assumed to be elastically restrained 
from rotations by restoring moments which are of equal strength on opposite edges and 
are proportional to the edge angle of rotation. 


Differential Equation and Boundary Conditions 

The small-deflection equilibrium equation for motion of an orthotropic panel in the 
presence of inplane tensile or compressive loads and supersonic flow may be written in 
the following form (ref. 16): 

D x 3 4 w + 2 L f M y D x \ 9 4 w [ °y 9 4 w 

1 ~ Px^y ^x 4 \ ^ 1 “ Px^y/ 9x 2 9y 2 1 - PxPy 9y 4 


+ N-v 


9 2 W 

ax 2 


+ N,; 


9y 


+ = z(x,y,t) 

9t Z 


(i) 


The panel bending stiffnesses in the x- and y-directions are denoted by D x and D y , 
respectively; D xy is the panel twisting stiffness, and y is the panel mass per unit 
area. The lateral loading induced by the supersonic flow is assumed to be given by two- 
dimensional quasi-steady aerodynamic theory so that (ref. 17) 

« = < 2 > 

where q = A-pU 2 is the free-stream dynamic pressure and /3 = \/m 2 - 1 is the 
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compi'essibility factor. The first term on the right-hand side of equation (2) corresponds 
to a static loading and the second term is a viscous loading corresponding to the aerody- 
namic damping. Reference 5 has shown that for M > 1.6 and from 0 to 10, the use 

D 

of two-dimensional static aerodynamic theory (aerodynamic damping neglected) yields 
flutter results in good agreement with those predicted by more exact aerodynamic theory. 


For finite rotational restraint and nondeflecting supports the boundary conditions 
(ref. 16) are as follows: 


Px 9^w 
1 ~ M-x^y 9x2 


9w 

9x 


- 9 X - — = 0 and w = 0 at x = - § 


1 ■ Mx^y 9x 2 


3 2 w + 0 |w = 0 ^ w = 0 at x=+ a 
x 9x 2 


Py 9 2 w e 9w 

1 - n xi x y 9y2 y 9y 

p y 9 2 w , c 9W 

1 - n x n y 3y2 y 9y 


0 and w = 0 at y = - ^ 
J 2 

> 

0 and w = 0 at y = 


(3a) 


(3b) 


where 6 X and Qy are the spring constants of the rotational springs on the boundaries 
x = and y = ±p respectively. 

The frequency-independent, linear, hysteretic formulation of structural damping is 
employed in the present analysis. Such a formulation may be derived for a one- 
dimensional case as follows . Assume that the material behaves as that for a Kelvin- 
Voigt viscoelastic body. For such a material the stress a is a function of strain rate 

— as well as of strain e: thus, 

8t 

a = kie + k 2 (4) 

ot 


where kj represents the elastic modulus E and k 2 is a coefficient of viscosity; both 
kj and k 2 are positive real quantities. 

For simple harmonic motion the strain may be written 


iwt 

e = e 0 e 

where £q is the maximum strain amplitude and oj is real. 


(5) 
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The stress now becomes 

°’= k l( 1 + iw ^) e 0 e 1Wt ( 6 ) 

ko 

where the term iu> — represents a stress 90° out of phase with the strain. For 
k l wko 

frequency-independent, linear, hysteretic structural damping, is replaced by a 

k l 

structural damping coefficient g which is independent of the magnitude of the frequency 
co. The stress may now be written 


a = E(1 + ig)e 0 e lwt (7) 

Thus, inclusion of this type of structural damping can be accomplished in dynamic analy- 
ses by using (1 + ig) to modify all terms which contain the elastic modulus E. This pro- 
cedure is most nearly correct for simple harmonic motion. Therefore, although this pro- 
cedure may not be strictly valid for all regions of the solution to the present problem, it 
should yield good results at the point of sustained simple harmonic motion which is 
assumed to be the point of incipient flutter for the panel. With this procedure, only the 
terms in equation (1) involving panel bending will be modified by a structural damping 
coefficient. In reference 15 plate vibration equations utilizing stress -strain relations 
similar to equation (4) also indicate that only those terms of the equations associated with 
bending should be modified by a structural damping coefficient. However, in the present 
analysis a structural damping coefficient associated with membrane force terms is also 
retained so that the present results can be compared with those of previous investigations. 
The damping is introduced by multiplying the bending and membrane force terms in equa- 
tion (1) by (l + igg) and (l + igj^, respectively, where gg and g M are the bending and 
membrane structural damping coefficients. Such a procedure is equivalent to assuming 
complex bending stiffnesses and for g M * 0, complex stress resultants. 


Solution 


Following the introduction of the generalized aerodynamic loading and the structural 

tt 4d x 


damping and dividing by 


1 " Mx^y 


, equation (1) becomes 


/I . \f 1 9 4 w . 2 Dl2 9 4 w , 1 d 2 9 4 w\ \ / k x 9 2 w , k y aV 

( gB V ax 4 7T 4 Dj 9 x 29 y 2 ,4 D x dy *J + ( 1§M )^2 b 2 9x 2 ^2 9y2 


4 . 1 92 w g A 9w A 9w _ q 

u> r 2 a 4 at 2 o) r a 4 8t 7T 4 a 3 8x 
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where 


D 1 


Px 

1 - Px/^y 


d 12 = D xy + 


(i,yD x 
1 “ PxPy 


(9a) 

(9b) 


D 2 = 


1 - 


^y 


k x = 


N x b 2 

7t2Di 


k y = 


0J r 


2 = 


Nyb^ 

7r2D^ 

7T^D^ 


yap 


(9c) 

(10a) 

(10b) 

(ID 


SA 


pc 

yco r 


( 12 ) 


2qa^ 

^r 


(13) 


A solution to equation (8) for sustained motion of the panel may be obtained by appli- 
cation of the Galerkin technique. The lateral deflection is a function of the spatial coordi- 
nates x and y and time t and may be expressed in the form 


w(x,y,t) =^A rg 


r s 



(14) 


where a is of the form 


a = cp + ico 


(15) 


The deflection functions Fr(a) anc * ^ s (b) 3X6 * aken to tke zero-stress vibration- 
mode shapes for supported uniform beams with equal elastic end rotational restraints and 
hence satisfy all the boundary conditions (eqs. (3a) and (3b)). On the basis of past experi- 
ence for a clamped panel (ref. 3), use of the beam vibration modes should yield excellent 
results. 
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Hedgepeth (ref. 18) has shown that for simple supports there is no stiffness coupling 
between cross-stream modes. In reference 19 it was indicated that such a coupling was 
insignificant for a fully clamped panel; therefore, coupling between cross-stream modes 
is assumed to be small for elastic restraints and is neglected herein. The two- 
dimensional aerodynamic theory precludes any aerodynamic coupling of orthogonal cross- 
stream modes. Thus, a one-term expansion in the cross-stream direction corresponding 
to the lowest beam vibration mode is used to determine the panel dynamic stability crite- 
rion. Substituting the assumed deflection (eq. (14)) into equation (8), multiplying by 

integrating over the area, and noting that for a nontrivial solution to equa- 
tion (8) theaeterminant of the coefficients of A rl must equal zero yields the following 
stability criterion: 


V 


' l + lg ^(?)0 +Bl ® +B2 
v*/\ o/n x J, n 


+ ( l + * b 4 


” 4 ( J 3 ) u 


A 31 

» 4 ( K 3) 22 


f + t^Bi + (l + ig M )B^| 


( l + *bb) 


7nAEi\ + B fe\ 

\A k 3 l w.. 


+ ( l + ig M)f 3 @ 22 + 


*4 


+ 4 / 


A L 32 
” 4 ( J *) M 


(l + ig B ) B l + (1 + ig M )B 3 - 


N, 


X L 23 

" 4 ( R 3) 22 


? ( J 3 ), 


|(l + ig B )B 1 + ( 1 + i g M )B3]^ 


( i + ib #® 33 + Bi ( 4 + H 


^ ( J 3 )„ 


7T 4 (k 3 ; 


-1 ( Kz ) 


" 4 M,. 


( 1 + ^ Pl ) 44 +Bi ®) 44 + B21 

+ ( i+ ^ M )te) 


zr) +B 4 
'44 J 


(16) 
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In the determinant of the coefficients of A r j, the coefficients (J i) mr > (Ki) m r > Vr’ 
and L mr represent integral expressions involving the streamwise deflection function 
and the constants Cq, C 4 , and C 2 represent integral expressions involving the cross 
stream deflection function. Details of the procedure and the integral expressions are 
given in the appendix. The coefficients Bi to B 4 are given as 

B = 2P 12 /a\ 2 C 1 
^ tt 2 Dj V 3 / TT^Cg 


b 2 


_ ^2 /a \ 4 C 2 


Dl V 3 / 7 T^Cq 


R _ k x /a\ 
Bs - 

Ba = k Jf 


J_c 1 
y W 


0 


(17) 


and 


\js = = 


a* 


00 r 


g A 


a 


(18) 


The expressions for B^ and B 2 contain the panel geometric and stiffness parameters 
and the expressions for B 3 and B 4 contain the midplane loading terms k x and k y . 
The coefficients Bj_, B 2 , and B 4 are also functions of the cross-stream boundary con- 
ditions. (See the appendix.) The term \p in equation (16) assumes the role of a complex 
eigenvalue which must be examined to determine the stability of the panel and includes 
both the panel frequencies and the aerodynamic damping coefficient. 

The panel behavior is characterized by the variation of cp + ico with an increase 
in aerodynamic load X. Instability of the panel occurs when cp becomes positive and, 
hence, the condition when cp = 0 is taken to be the critical condition. This condition 
corresponds to sustained simple harmonic motion and represents incipient flutter for the 
panel. The relations necessary to determine cp and 00 may be obtained by solving for 
a from equation (18), as follows: 



T^l/X 2 * 4 (*r + %) 


(19) 


Equation (19) yields two roots for a and, thus, two sets of (cp + iw). The roots for co 
occur as ± the same absolute quantity; however, examination of equations ( 6 ) and (7) 
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reveals that introduction of positive structural damping in the flutter equation implies that 
the algebraic sign of the damping coefficient g and oo be the same. Thus, for a posi- 
tive damping coefficient g the solution is valid only for the root of a corresponding to 
positive values of co. For certain panel conditions if both roots for a are considered, 
a false instability (c p > 0) is found to occur over a wide range of X including X = 0, a 
completely unrealistic result. The anomalies encountered in considering the negative 
root for co have been noted by other authors in references 20 and 21. 

Substituting equation (15) for a into equation (19), equating the real and imaginary 
parts, and selecting the root corresponding to positive co yields 

<P _ g A , n /P + gA 2 + 4 ^R 
co r 2 2 j 2^j 2 

and 

co _1 /P ~ g A 2 ~ 4 ^R 
C0 r 2 If 2 

where 

P = j/(g A 2 + 4* r ) 2 + (4 

Equations (20) and (21) govern both the dynamic and static stability of the panel. Zero 
values of the term under the radical in equation (21) yield zero values of the frequency 
which indicate zero stiffness and hence loss of static stability for the panel. Since cp = 0 
is sought as the critical condition for dynamic instability, equation (20) is used to deter- 
mine the relation between g^ and the real and imaginary parts of the eigenvalue at 
flutter. The aerodynamic damping coefficient is defined as a positive quantity; 
therefore, only positive values of i//j in equation (20) can lead to dynamic instability. 

The relation between g^ and the real and imaginary parts of the eigenvalue thus 
obtained is 



As both the aerodynamic damping and the structural damping approach zero, the use 
of equation (23) leads to values of X which approach the values required for coalescence 
of the two lowest natural frequencies. In the presence of damping, flutter occurs at a 
value of X corresponding to the value required to cause cp to become positive. For 
specific flow conditions and panel material, the aerodynamic damping coefficient may be 
expressed as a function of X by combining equations (11) to (13) and is given by 


( 20 ) 


( 21 ) 


( 22 ) 
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( 24 ) 





Thus, to obtain flutter boundaries a solution for the eigenvalue <|/ as a function of 
the aerodynamic loading parameter A was programed for a digital computer and utilized 
an eigenvalue routine for a square complex matrix. To obtain points of the flutter bound- 
ary for a particular panel, it is necessary to specify the panel geometric, stiffness, mid- 
plane loading, and structural damping characteristics as well as the degree of elastic 
rotational restraint present at the panel boundaries. The aerodynamic loading param- 
eter X is then varied until both equations (23) and (24) yield identical values of g A . 
Limitations of the computer program restricted the maximum number of streamwise 
terms in the assumed deflection (eq. (14)) to 24. Hence, the program yields accurate 
results only for problems in which 24 terms are sufficient for convergence. 


RESULTS AND DISCUSSION 


Data are presented to evaluate the accuracy of the effects of aerodynamic damping 
given herein by the two-dimensional quasi-steady approximate aerodynamic theory as 
compared with the effects given in reference 5 by the more exact three-dimensional 
unsteady aerodynamic theory, to establish the validity of associating linear hysteretic 
structural damping with panel bending only, and to determine whether such a procedure 
yields a better correlation between theoretical and experimental panel flutter results. 
Although the analysis was developed in a general form to include orthotropic panels, 
results are presented for isotropic panels only. Comparison of representative results 
from the present analysis with results obtained from the exact analysis presented in 
reference 1 revealed that the present results were accurate for clamped panels having a 
length-width ratio less than about 8. 

Damping Obtained From Aerodynamic Theories 

Flutter results for a simply supported panel obtained from both two-dimensional 
quasi-steady aerodynamic theory and three-dimensional linearized supersonic potential 
flow theory presented in reference 14 and utilized in reference 5 are given in figure 2 for 
a panel with a length-width ratio of 4. The flutter parameter A^/^ is shown as a func- 
tion of k X/ /k x cr , the ratio of the midplane compressive load to the critical value required 
for buckling with no airflow. The dash-line curve represents the flutter boundary obtained 
from two-dimensional static aerodynamic theory (no damping). The numbers on this curve 
indicate the modes that coalesced for flutter. The circles represent flutter results from 
the three-dimensional potential flow theory taken from figure 9 of reference 5, which were 
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calculated from a six mode solution for aluminum panels at sea level for M = 3.0; this 
value of Mach number corresponds to the value at which most published experimental data 
for stressed panels have been obtained. The solid-line curve represents results from 
two-dimensional aerodynamic theory and was obtained by utilizing 10 modes in the pres- 
ent analysis. Structural damping was zero for all cases and for the present analysis aero- 
dynamic damping is given by equation (24). The fact that the buckling load in the presence 

of airflow can be larger than that for no airflow causes the curves to extend beyond 
k x 

— — = 1.0. The buckling point on the flutter boundary in the presence of airflow is 
^x,cr 

defined as the transition point. For no damping, the boundary indicates that zero values 
of the flutter parameter occur at values of k X/ A x?cr of approximately 0.58, 0.71, and 
0.89. When the aerodynamic damping g^ is included in the calculations, the effect is 
the removal of the zero-dynamic-pressure points although the saw-toothed-like character 
of the boundary remains. As can be seen, the differences between the results from the 
two-dimensional quasi-steady aerodynamic theory and the three-dimensional potential 
flow theory are slight and may be considered insignificant in view of the fact that the lat- 
ter results were obtained by a method which entailed a great deal of cross plotting of 
results from a six mode solution, which were probably not as well converged as the pres- 
ent results. Reference 5 shows that results from static aerodynamic theory agree rea- 
sonably well with those from three-dimensional potential flow theory for Mach numbers 
greater than 1.6; therefore, it is assumed that results from two-dimensional quasi-steady 
aerodynamic theory will also agree for M > 1.6. Two-dimensional quasi-steady aerody- 
namic theory does yield satisfactory results, when compared with results from three- 
dimensional potential flow theory, and provides a greatly simplified solution for many 
aeroelastic problems in which the effects of aerodynamic damping are important. 

Effects of Structural Damping 

Figure 3 illustrates the effects of structural damping on the flutter boundary pre- 
sented in figure 2. The curves represent results from the present analysis for values of 
the aerodynamic damping coefficient from equation (24) for aluminum -alloy panels at sea 
level for M = 3.0. The dot-dash-line curve for no structural damping is taken from fig- 
ure 2. The dash- line curve was obtained for a value of the bending and membrane struc- 
tural damping coefficients equal to 0.01. The solid-line curve represents results for a 
value of the bending damping coefficient of 0.01 and a value of the membrane damping 
coefficient of zero. As shown by the dash-line curve, structural damping with equal 
bending and membrane damping coefficients tends to smooth out the saw-toothed-like 
flutter boundary; however, this effect lessens as k X/ /k xcr approaches the transition 
point. Similar results were obtained in reference 5 which used equal structural damping 
coefficients. Thus, this type of structural damping appears to have little effect on the 
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flutter of panels stressed to the point of buckling (transition point). However, as shown by 
the solid- line curve, when membrane damping is taken to be zero, structural damping 
smooths out the boundary to a greater degree and also raises the boundary at the point of 
buckling. 

Since the effects of aerodynamic damping and structural damping with g-g = g-^ 
appear to be small as k x /k x ,cr approaches the transition point (see fig. 3), it is of inter- 
est to examine a special case for which the theory, damping neglected, predicts flutter at 
zero dynamic pressure at the panel buckling load. Such a condition is shown in figure 4 

wherein the flutter parameter is plotted as a function of k x /k x cr for a rota- 

' N v 

tionally restrained panel with a length-width ratio of 3.3 and a stress ratio of — = 1. 

N x 

For this condition k x>cr is equivalent to k^-p, the transition-point value of k x . The 
lower solid-line curve represents the flutter boundary for no damping, and the dash-line 
curve represents the flutter boundary for gg = gj^ = 0.01 and values of g^ obtained 
from equation (24) for an assumed aluminum -alloy panel at Mach 3 and sea-level flow 
conditions. The upper solid-line curve is the flutter boundary obtained for the same 
flow conditions but with gj^ = 0. As was shown in figure 3, if gg = gj^, neither struc- 
tural damping nor aerodynamic damping has any effect at the transition point (buckling 
point). Thus, for this special case (dash-line curve), theory predicts the occurrence of 
flutter at the transition point for zero dynamic pressure. This anomaly completely dis- 
appears when structural damping is associated with bending only. The trend of the flutter 
boundary thus obtained is physically reasonable and is typical of existing experimental 
boundaries. Therefore, it appears that structural damping is more realistically repre- 
sented by modifying only the bending terms of the differential equation through use of 
complex stiffness coefficients. 

Effects of Edge Rotational Restraint, Stress Ratio, and 
Length-Width Ratio on Transition Point Flutter 

The structural damping results are particularly significant when it is realized that 
existing experimental results indicate that the most critical flutter condition occurs near 
or at the buckling point of the panel. This condition as defined by experiment is the tran- 
sition point (intersection of prebuckled panel boundary and postbuckled panel boundary) 
and, as indicated in figure 4, past attempts at predicting the minimum value have in many 
instances led to the physically unreasonable zero values of dynamic pressure. Addition- 
ally, theoretical results which neglect damping or include structural damping in the form 
gg = gjyi indicate that the transition point (k xcr = k x> T) is extremely sensitive to varia- 
tions in edge rotational restraint, stress ratio, and length-width ratio. Therefore, results 
using the present representation of structural damping are presented concerning the effect 
of these parameters on transition point flutter. 
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Figure 5 illustrates the effect of varying edge rotational restraint for a panel having 
a length-width ratio of 3.3, a stress ratio of 1, and equal rotational springs on all edges. 
The flutter parameter at the transition point A-p 1 / 3 is plotted as a function of the edge 
rotational restraint coefficient q x . The lower curve represents results for any value of 
aerodynamic damping and any equal values of bending and membrane structural damping 
since dam ping included in this manner has no effect on flutter. The upper curve repre- 
sents results for any value of aerodynamic damping and a value of bending structural 
damping of 0.01. A decrease in flutter resistance with an increase in restraint is exhib- 
ited by both curves. (This result agrees with ref. 1 wherein it was found that, for no 
damping, as edge rotational restraint increases, a stressed panel can become more sus- 
ceptible to flutter even though its resistance to buckling increases.) In fact, the lower 
curve in figure 5 exhibits, at a value of q x = 40, a zero-dynamic-pressure transition 
point. Not only does the present analysis eliminate this anomaly, as shown by the upper 
curve, it also predicts a less severe decrease in flutter margin with an increase in edge 
rotational restraint. 

The effects of varying stress ratio on transition point flutter are illustrated in fig- 
ure 6 for a panel having a length-width ratio of 4 and clamped edges. Again, the lower 
curve represents results which are unaffected by aerodynamic damping and structural 
damping when gg = gjyj, and the upper curve represents results obtained by associating 
structural damping with bending only. Comparison of the two curves reveals that asso- 
ciating structural damping with bending only not only removes the many zero-dynamic- 
pressure flutter points (which occur when the panel has an equal choice of buckling modes) 
but also gives a nearly constant value of the transition point flutter parameter over a wide 
range of stress ratio. 

Figure 7 illustrates the effects of varying length-width ratio on the transition point 
flutter for a panel with clamped edges and a stress ratio Ny/N x of 1. The lower curve 
represents results unaffected by damping, and the upper curve is again representative of 
results associated with bending structural damping only. Comparison of the two curves 
reveals that associating structural damping with bending only removes the extreme sensi- 
tivity of the transition point to variations in length-width ratio. 

The results of figures 5, 6, and 7 suggest that the apparent sensitivity of critical 
flutter boundaries to edge rotational restraint, stress ratio, and length-width ratio was 
spurious and was associated largely with the anomalous zero values of critical dynamic 
pressure which, in turn, were due to an unrealistic theoretical model of structural 
damping. 
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Comparison of Theory and Experiment 

The experimental flutter boundary presented in reference 2 for a panel with a length- 
width ratio of 2.9 and an average edge rotational restraint coefficient q x of 44 is com- 
pared with theoretical boundaries obtained from the present analysis in figure 8. Only 
the portion corresponding to stresses below buckling is shown. Beyond buckling, the 
experimental boundary rises and small deflection theory is not sufficient to handle this 
condition. The boundaries are presented in terms of the flutter parameter and 

the percent of midplane compressive load required for buckling in the presence of airflow. 
Theoretical boundaries are represented by the solid- line curves; the lower curve is for 
zero damping and the upper curve is based on an estimated bending damping coefficient 
gg of 0.01 and values of the aerodynamic damping coefficient g^ calculated from the 
experimental test conditions. Reference 22 indicates that the maximum value of material 
damping for an aluminum alloy similar to that used in reference 2 is approximately 0.003. 
References 23 and 24 reveal that damping mechanisms present at panel boundaries may 
increase the structural damping up to five times the value for material damping. It is, 
therefore, believed that the estimated value of the bending damping coefficient of 0.01 is a 
realistic value for the experimental data represented by the circular symbols. 

Agreement between the experimental boundary and the theoretical boundary for no 
damping is reasonably good for low values of midplane compressive stress but becomes 
poor in the region of the transition point. Aerodynamic damping and structural damping 
associated with bending only improves the agreement at the transition point, but the theory 
remains conservative. The values of q x for the points shown range from about 27 to °° 
or fully clamped. Therefore, it is believed that the experimental boundary is not truly 
representative of a boundary for q x = 44 and that some scatter should be expected. 

Thus, a need exists for experimental flutter boundaries wherein the variation of q x is 
not so large to provide data for a more adequate evaluation of the present theory. 

CONCLUSIONS 

The flutter of stressed panels with elastic edge rotational restraint was investigated 
theoretically. A modal solution which includes both structural and aerodynamic damping 
was presented for the supersonic flutter of flat orthotropic panels. The solution utilizes 
two-dimensional quasi-steady aerodynamics and is valid for panels with nondeflecting 
edges . The theoretical results and the comparison of a theoretical boundary from the 
present analysis with an experimental boundary presented in NASA TN D-3498 revealed 
that the following conclusions can be made : 
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1. Two-dimensional quasi-steady aerodynamic theory yields satisfactory results 
when compared with three-dimensional potential flow theory for M > 1.6 and may pro- 
vide greatly simplified solutions to aeroelastic problems wherein the effects of aerody- 
namic damping are important. 

2. Inclusion of structural damping in a small deflection plate theory is more realis- 
tically accomplished by modifying only the bending terms of the differential equation 
through the use of complex bending stiffness coefficients; such a representation removes 
all physically untenable results that otherwise occur and provides reasonable agreement 
between theoretical and experimental panel flutter results over the entire stress range 
from zero to buckling. 

3. Representation of structural damping by use of complex bending stiffness coeffi- 
cients renders the flutter parameter at the transition point (transition from flat to buckled 
panel) relatively insensitive to variations in edge rotational restraint, stress ratio, and 
length-width ratio. 

The ability to predict experimental results with reasonable accuracy over the entire 
stress range makes it possible to consider the flutter design of stressed panels on a ratio- 
nal analytical basis rather than continuing to rely on empirical boundaries which must be 
periodically revised as new experimental data are obtained. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., September 23, 1968, 

126-14-02-22-23. 
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APPENDIX 


SOLUTION TO FLUTTER EQUATION 


A solution to the flutter equation (eq. (8)) may be obtained by application of the 
Galerkin technique. Substitution of the assumed deflection given by equation (14) into 

equation (8), multiplying by F m(a) G l(b)’ and inte S ratin g over the area results in the 
following set of equations for the amplitude coefficients A r j: 
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The primes and superscript Roman numeral denote differentiation with respect to — or 
y a 

b’ 

The integrals in equations (A3) to (A 6 ) are evaluated in terms of the expressions on 
the right-hand sides of these equations by use of the following deflection functions for the 
zero stress vibrations of beams with equal end rotational restraint (see ref. 25): 


Symmetrical modes 



__ cos 61 

cos 2a i 5 - cosh 26 i I 


(J = 1, 3, 5 . . .) (A7) 


Asymmetrical modes 

/■y \ „ sin 6 i 

Fj(a) = sin 2 5j | - sinh 2 6j | (j = 2, 4, 6 . . .) (A8) 

where is a constant dependent on the degree of rotational restraint. The expressions 
for determining 6 j from reference 25 in the terminology of the present analysis are as 
follows : 
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Symmetrical modes 
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Asymmetrical modes 
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(All) 


The integral expressions in equations (A 2 ) may be evaluated by use of equations (A7) to 
(A9), that is, replacing ^ with replacing q x with q y = ^2., and taking j = 1. 
Expressions for the integrals in equations (A3) to (A 6 ) are given as follows: 
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For a nontrivial solution to equation (Al) the determinant of the coefficients of A r i 
must be zero. The determinant is given in the main text (eq. (16)). 
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